Random antiferromagnetic quantum spin chains: 
Exact results from scaling of rare regions 



Ferenc Igloi^^^, Robert Juhasz^'^ and Heiko Rieger'^ 
^ Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 
^ Institute for Theoretical Physics, Szeged University, 11-6720 Szeged, Hungary 
^ FB 10.1 Theoretische Physik, Universitdt des Saarlandes, 6604-1 Saarhriicken, Germany 

(April 11, 1999) 

We study XY and dimerized XX spin-1/2 chains with random exchange couphngs by analytical 
and numerical methods and scaling considerations. We extend previous investigations to dynamical 
properties, to surface quantities and operator profiles, and give a detailed analysis of the Griffiths 
phase. We present a phenomenological scaling theory of average quantities based on the scaling 
properties of rare regions, in which the distribution of the couplings follows a surviving random 
walk character. Using this theory we have obtained the complete set of critical decay exponents of 
the random XY and XX models, both in the volume and at the surface. The scaling results are 
confronted with numerical calculations based on a mapping to free fermions, which then lead to an 
exact correspondence with directed walks. The numerically calculated critical operator profiles on 
large finite systems (L < 512) are found to follow conformal predictions with the decay exponents 
of the phenomenological scaling theory. Dynamical correlations in the critical state are in average 
logarithmically slow and their distribution show multi-scaling character. In the Griffiths phase, 
which is an extended part of the off-critical region average autocorrelations have a power-law form 
with a non-universal decay exponent, which is analytically calculated. We note on extensions of our 
work to the random antiferromagnetic XXZ chain and to higher dimensions. 



I. INTRODUCTION 

Quantum spin chains exhibit many interesting phys- 
ical properties at low temperatures which are related 
to the behavior of their ground state and low-lying ex- 
citations. In this context one should mention quasi- 
long-range-order (QLRO), topological order and quan- 
tum phase transitions, which have purely quantum me- 
chanical origin. Considering isotropic antiferromagnetic 
chains for integer spins there is a gap, whereas half inte- 
ger spin chains are gap-less |^|. However, alternating cou- 
plings in spin-1/2 chains yield a dimerized ground state 
that has physical properties similar to the spin-1 chain: 
there is a finite gap, spatial correlations decay exponen- 
tially and there is string topological order. 

Randomness may have a profound effect on the phys- 
ical properties of quantum spin chains, as demonstrated 
by recent analytical and numerical studies ||^. As an 
interplay of randomness and quantum fluctuations there 
are new exotic phases in disordered quantum spin chains, 
which are not present in classical random or pure quan- 
tum systems. It has been noticed, that pure gap-less 
systems are generally unstable against weak randomness 
whereas for gaped systems a finite amount of dis- 
order is necessary to destroy the gap (but see also 

i). 

Among the theoretical methods developed for disor- 
dered quantum spin chains one very powerful procedure 
is the renormalization group (RG) approach introduced 
by Dasgupta and Ma ||^. This RG method, which is 
expected to be asymptotically exact at large scales, i.e. 



close to critical points, has been applied for a number 
of random quantum systems. The fixed point distribu- 
tion of the RG transformation has been obtained ana- 
lytically for some random quantum spin chains, among 
others for the transverse Ising spin chain the spin- 
1/2 Heisenberg and related spin chains with random an- 
tiferromagnetic couplings 0. On the other hand some 
other one-dimensional problems {S = 1/2 Heisenberg 
chain with mixed ferromagnetic and antiferromagnetic 
couplings ||l^,|ri|, S = 1 antiferromagnetic chain with 
and without biquadratic exchange etc.), as well as 

higher dimensional random quantum systems |T2| have 
been studied by numerical implementation of the RG pro- 
cedure. Comparing the RG results with those obtained 
by direct numerical evaluation of the singular quantities 
p3|-p^ and by other exact and numerical meth- 

ods one has obtained a good agreement in the vicinity of 
the critical point. 

There are, however, other interesting singular quanti- 
ties, which are not accessible by the RG method. We 
mention, among others, the dynamical correlations |l^ ] 
and the behavior of the system far away from the crit- 
ical point in the Griffiths phase jl^], which denotes an 
extended region of the parameter space around the crit- 
ical point. In the Griffiths phase the system is gap- 
less, thus dynamical correlations decay with a power-law, 
however there is long-range-order with exponentially de- 
caying spatial correlations. For the random quantum 
Ising chain dynamical correlations, both at the critical 
point and in the Griffiths phase have been exactly deter- 
mined pG E2| using a mapping to the Sinai model p3[. 
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i.e. random walk in a random environment. 

In this paper we are going to study the S = \/2 dis- 
ordered XX and XY spin chains by analytical and nu- 
merical methods and by phenomenological scaling theory. 
The RG treatment of the problem by Fisher predicts 
the antiferromagnetic random XX fixed point to control 
the critical behavior of the antiferromagnetic Heisenberg 
{XXX) model, too. Furthermore, for random isotropic 
chains the RG approach predicts QLRO, thus the av- 
erage spatial correlations of different components of the 
spin decay with a power-law. In this so called random 
singlet (RS) phase all spins are paired and form singlets, 
however, the distance between the two spins in a singlet 
pair can be arbitrarily large. Then these weakly cou- 
pled singlets dominate the average correlation function, 
therefore all components of the correlation function are 
predicted to decay with the same exponent. 

Leaving the critical state by introducing either 
anisotropy or dimerization, randomness will drive the 
system into the Griffiths phase. As shown by an RG anal- 
ysis 1^, applicable in the vicinity of the RS fixed point, 
the Griffiths phase is characterized by the dynamical ex- 
ponent, 0, defined by the asymptotic relation between 
relevant time {tj.) and length scales (^) as 



(1.1) 



The dynamical exponent is predicted to be a continuous 
function of the quantum control parameter (anisotropy 
or dimerization) and the singular behavior of different 
physical quantities (specific heat, susceptibility, etc.) are 
all expected to be related to the value of the dynamical 
exponent. 

The RG predictions by Fisher Q and others have 
been confronted with results of numerical studies p4]-p6[ , 
especially in the RS phase of isotropic chains, but some 
cross-over functions of correlations have also been studied 
in the Griffiths phase. In the RS phase some numerical 
results are controversial: in earlier studies psf a different 
scenario from the RG picture is proposed (in particular 
with respect to the transverse correlation function) , later 
investigations on larger finite systems have found satis- 
factory agreement with the RG predictions |^ , although 
the finite-size effects were still very strong. 

In the present paper we extend previous work in several 
directions. Here we consider open chains and study both 
bulk and surface quantities, as well as end-to-end corre- 
lations. We develop a phenomenological theory which is 
based on the scaling properties of rare events and deter- 
mine the complete set of critical decay exponents. We 
calculate numerically (off-diagonal) spin-operator pro- 
files, whose scaling properties are related to (bulk and 
surface) decay exponents [ p7| and compare the profiles 
with predictions of conformal invariance. Another new 
feature of our work is the study of dynamical correla- 
tions, both at the critical point and in the Griffiths phase. 
Finally, we perform a detailed analytical and numerical 
study of the Griffiths phase and calculate, among others, 
the exact value of the dynamical exponent, z in Eq.(l.l). 



The structure of the paper is the following. The model 
and its free-fermion representation are presented in Sec- 
tion 2. A phenomenological theory based on the scaling 
behavior of rare events is developed in Section 3. Re- 
sults in the critical state, where there is quasi-long-range 
order in the chains is presented in Section 4, whereas 
the Griffiths phase is studied in Section 5. We discuss 
the extensions of our results to random antiferromagnetic 
XXZ chains and to higher dimensions in the final Sec- 
tion, whereas some technical calculations are presented 
in the Appendices. 



II. THE MODEL AND ITS FREE-FERMION 
REPRESENTATION 

A. The XY and XX models 

We consider an open XY chain (i.e. with free bound- 
ary conditions) with L sites described by the Hamilto- 



L-l 

1=1 



TV cv qv ^ 



(2.1) 



where the Sj^ {ii — x,y) are spin-1/2 operators and the 
couplings {Jj^ > 0) are independent random variables 
with distributions 7r'^(J^). The quantum control param- 
eter is the average anisotropy defined as: 



Sa 



[In 



[InJ^la 



var[ln J^] + var[ln Jv] 



(2.2) 



where var(a;) is the variance of random variable x and 
[. . .]av denotes average over quenched disorder. For 6a > 
(< 0) there is long-range-order in the x (y) direction, 
i.e. limr^oo[C^(7')]av 7^ 0, where 



(2.3) 



and for 6a ~ the system is in a critical state with quasi- 
long-range-order, where correlations decay algebraically, 
i.e. 



(2.4) 



In the XX model, where the x and y couplings are cor- 
related as = Jf = Ji, we introduce alternation such 
that even (e) and odd (o) couplings, connecting the site 
2i, 2i -I- 1 and 2i — 1, 2i, respectively, are taken from dis- 
tributions p'^(Je) and p°{Jo), respectively. For the XX 
model the quantum control parameter is the average 
dimerization defined as: 



6d 



[In Jo] av ~ [In Je] av 
var[ln Jo] + var[ln Jg 



(2.5) 



The RS phase is at i5d = 0, whereas 6d corresponds 
to the random dimer (RD) phase. Throughout the paper 
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we use two types of random distributions, both for the 
XY and XX models. For the XY model with the binary 
distribution the couplings can take two values A > 1 
and 1/A with probability p and q — 1 — p, respectively, 
while the couplings are constant: 

7r^(J^) = pS{.r - A) + q6{r - x-^' 



TTy{jy) = 5{jy - 4) 



(2.6) 



At the critical point [p — q)\n\ = In Jq. The uniform 
distribution is defined via 



1, for < < 1 
0, otherwise 

0, otherwise 



for o<jy <jy 



(2.7) 



and the critical point is at Jq = 1. 

For the XX model the corresponding distributions 
p'^{J'^) and p°{J°) follows from the correspondences: 



r J^ jy j°, 
7r^(j^) ^ p^iJ"), TTy{jy) p°{j°) 



(2., 



Note that the critical points of the two models {5a — 
and 5d = 0, respectively) are not equivalent due to the 
different disorder correlations. 



B. The XY chain and the directed walk model 



quadratic form in fermion operators. It is then diago- 
nalized through a canonical transformation which gives 



L 1 



(2.9) 



The fermion excitations are non-negative and satisfy the 
set of equations 



e,*,(0 = Jf_i$,(;-1) + Jr'i>,(; + 1) 
6,ci>,(0 = Jr_i*,(/-1) + Jf*,(? + 1); 



(2.10) 



with the boundary conditions Jf, = — 0. The vec- 
tors $q's and 5'g's which are related to the coefficients of 
the canonical transformation are normalized. They enter 
into the expressions of correlation functions and thermo- 
dynamic quantities. 

Usually one p roceeds p8| by eliminating either ^'^ or 



$q in Eqs.(2.1C) and the excitations are deduced from 



the solution of quadratic equations. This last step can 
be avoided by introducing a 2L-dimensional vector Vq 
with components: 

Vq{Al - 3) = $,(2? - 1), Vq{Al - 2) - *,(2; - 1), 
Vq{Al-l) = ^q{2l), Vq{Al) = <^q{2l)- (2.11) 



Using the Jordan- Wign er t ransformation, the XY 
model Hamiltonian in Eq.(pT[) can be rewritten as a 



and noticing that the relations in Eqs.( 2T0| ) then corre- 
spond to the eigenvalue problem of the matrix: 
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(2.12) 



The matrix T can be interpreted as the transfer matrix 
(TM) of a directed walk (DW) problem on four inter- 
penetrating, diagonally layered square lattices. Each 
walker makes steps with weights Jf and Jf between next- 
neighbor sites on one of the four square lattices and the 
walk is directed in the diagonal direction (see Fig. |l|). 

According to Eqs.(2.1C), changing $q into in Vq, 



the eigenvector corresponding to — is obtained. Thus 
all information about the DW and the XY model is con- 
tained in that part of the spectrum with > 0. Later 



on we shall consider this sector. We note that similar 
correspondence has been established earlier between the 
DW and the transverse-field Ising model (TIM) ||§ . 

The eigenvalues of T in ( 2.12 ) are of two classes. For 
q = 2i — 1, i = 1, 2, . . . , L the odd components of the 
eigenvectors are zero, i.e. V2i-i(2j) = 0, j = 1, 2, . . . , L, 
whereas for the other class with q = 2i the even compo- 
nents are zero, V2i{2j — 1) = 0. Consequently T can be 
expressed as a direct product T — T^^Tr, where the 
trigonal matrices , Tr of size L x L represent transfer 
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matrices of directed walks. As a result one has to diago- 
nalize these two matrices of size L x L. Thus for chains 
with even number of sites, L = 2N, the two classes of 
eigenvectors are given in terms of the variables $ and ^' 
via: 

e2k-i- $2fc-i(2j) = 1-2^-1(2^-1) = , 
e2k. *2fe(2j-l) = *2fc(2j) = ^^■''^> 

for i, j = 1, . . . , TV. Furthermore we assume that the vec- 
tors $^and are normalized to 1 separately. 

For the XX model the even and odd sectors are degen- 
erate, e2fc_i — e2ki thus it is sufficient to diagonalize only 
one matrix. In this case one has the additional relations: 

$2fe-i(2j-l) - *2fc(2j-l) 
*2fe(2j) - *2fc-l(2j) . 



XX -model (2.14) 



The matrices and Tr are in one-to-one correspondence 
with the eigenvalue problem of one-dimensional TIM-s. 
This exact mapping for finite open chains is presented in 
Appendix A. 
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FIG. 1. Sketch of the directed walk prob lem correspond- 
ing to the transfer matrix given in eq. (2.12). Note that one 
has two independent walks in the t-direction of the diagonally 
layered square lattice, corresponding to the independent sub- 
spaces for the eigenvalue problem (see text). The coupling 
strength J^'^ are the transition rates for the random walker 
from one sites in row i to those in row i ± 1. 



C. Local order- parameters 

Next we are going to study the long-range-order in the 
ground state of the system. Having free boundary condi- 
tions, as in (2.1), the expectation value of the local spin 
operator (0|5f |0) (and (OjS'f |0)) is zero for finite chains. 
Then the scaling behavior of the spin operator can be ob- 
tained from the asymptotic behavior of the (imaginary) 
time-time correlation function: 



Gf (r) = {o\sr{T)srmo) 

("I 



|(n|^nO)Pexphr(i?„-£;o)] 



(2.15) 



where |0) and |n) denote th e gr ound state and the nth 
excited state of H in Eq. (2T), with energies Eq and 
En, respectively. In the phase with long-range-order the 
first excited state is asymptotically degenerate with the 
ground state in the thermodynamic limit, thus the sum 
in Eq. ( 2.15| ) is dominated by the first term. In the large 
T limit limT-^oo (t) = (mf)^, thus the local order- 
parameter is given by the off-diagonal matrix element: 



mf = {l\Sm . 



(2.16) 



In the free fermion representation S[ is expressed as l2S 



with 



(2.17) 



(2.18) 



Using |1) = rii\0) the matrix-element in Eq. ( 2.16| ) 
is evaluated by Wick's theorem. Since for i ^ 
j (0|A,Aj|0) = {0\BiBj\0) = we obtain for the local 
order-parameter : 



Hi Gil G12 
H2 G21 G22 



Hi Gil G 



12 



Gii^i 
G21-1 



Gii-i 



(2.19) 



where 



Gjk 



(0^1^,10) = 
mkA,\0) 



*i(j) 

.-^v|/^(fc)$,(j) 



(2.20) 



For surface spins the local order-parameter is simply 
given by mf = <I>i(l)/2, which can be evaluated in the 
thermodynamic limit L ^ 00 in the phase with long- 
range-order, when ei = 0. Using the normalization con- 
dition J2i I'J'iWP = 1 we obtain for the surface order 
parameter: 



E n 0^ 



L/2-1 I 

1+ E n 

1=1 j=i 



-1/2 



-1/2 



XY 



XX. 



(2.21) 



We note that this formula is exact for finite chains if 
we use fixed spin boundary condition, Sf^ = ±1/2, wich 
amounts to have J^_i = 0. In the fermionic description 
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the two-fold degeneracy of the energy levels, correspond- 
ing to S"! = 1/2 and Sf^ = —1/2, is manifested by a zero 
energy mode in Eq. (2.9) and from the corresponding 
eigenvector one obtains mf in Eq. ( 2.21 ) for any finite 
chain. 

For non-surface s pins the expression of the local order- 
parameter i n Eq. ( ^.19 ) can be simplified by using the 
relations in ( 2.13 ). Then, half of the elements of the de- 
terminant in ( 2.ig| ) are zero, the non-zero-elements being 
arranged in a checker-board pattern, and mf can be ex- 
pressed as a product of two determinants of half-size, 
which reads for I = 2j as: 



Hi 


Gl.2 


Gi_4 


Gl,2]-2 


Hs 


Gz^2 


GsA 


G3aj-2 


H2j-1 


G2j-1.,2 


G2j-1A 


■ ■ G2j-l,2j-2 


G2,l 


G2,3 


G2,2j-1 




G4,l 


G4,3 


Gi.2j-1 


(2.22) 


G2jA 


G2j.,3 ■ ■ 


G2j,2j-1 





The local order-parameter m]' , related to the off-diagonal 
matri x-elem ent o f the operator Sf can be obtained from 
Eqs. (|2.19D and ( |2.2lD by exchanging Jf 



For the S[ operator the autocorrelation function, 
Gf (t) , can be expressed in a similar way as Gfir) in Eq. 
( ^.15 ) and its long time Hmit, Hmi-^oo Gf{T) ~ {rnf)"^, is 
given by the local order parameter 



mf = (0,|5f|O) 



(2.23) 



Here \(j)z) denotes the lowest eigenstate of H in Eq. (2.1) 



having a non-vanishing matrix-element of S[ with the 
ground state. In the free fermion representation S[ can 
be written as (28| 

S[ = ^ABi (2.24) 
and the off-diagonal order-parameter is given by: 

mi=l\- $i(0*2(0 + *i(0*2(OI • (2.25) 



For the XX model one ca n ob tain simple expressions 
using the relations in Eqs. ( 2.14 ) as: 



^^,,_,^^[^l{2^-^)r 



XX~-mode\ 



(2.26) 



'■2i 



= mli2^)f 



D. Autocorrelations 

Next we consider the dynamical correlations of the sys- 
tem as a function of the imaginary time r. First, we note 
that the correlations between a;-component s of the sur- 
face spins can be obtained directly from eq(2.15) as: 



GUr) = lY.\^,il)\'eM-re,) 
q 

L/2 

= lEl*2«-i(l)l'exp(-Te2.-i), (2.27) 



where we have used the relations in Eq. ( ^.13| ). 

For bulk spins the matrix- element {n\Sf\Q) in eq( 2.15| ) 
is more complicat ed to evaluate, one has to go back to the 
first equation of (2.15) and considers the time-evolution 
in the Heisenberg picture: 

^r(T) = exp(rff)5fexp(-ri/) 

= ^Ai{t)Bi{t) . . . Ai^i{t)Bi^i{t)Mt) . (2.28) 

The general time and position dependent correlation 
function 

{Sr{r)Sr+„) = \{Ai{t)Bi{t) ■ ■ ■ A{t)AiBi . . . A,+„) , 

(2.29) 

can then be expanded using Wick's theorem into a sum 
over products of two-operator expectation values, which 
can be expressed in a compact form as a PfafRan: 



{Ai{r)Bi{r)) 



{Ai(r)A2{T)) 
{Bi(r)A2{T)) 



{Ai{r)B2{T)) 
{Bi{t)B2{t)) 
{Mr)B2iT)) 



{Ai{t)Ai{t)) {Ai{t)Ai) 
{Bi{t)Ai{t)) {Bi{t)Ai) 
{A2{t)Ai{t)) {A2{t)Ai) 



{Ai{T)Ai+n) 

{Bi{T)Ai+n) 
(A2(r)Ai+„) 



±[detCy]^/^ 



{Bl+n-lAl+n) 



(2.30) 
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where Cy is an antisymmetri c mat rix Cij = —Cji, with 
the elements of the Pfafhan ( 2.3C ) above the diagonal. 
At zero temperature the elements of the PfafRan are the 
following: 



{Aj{T)Ak) = exp(-Te, 
q 

{A,{T)Bk) = ^$,(j)*,(fc)exp(-Te,) , 
q 

{B,{T)Bk) = -^*^(j)*,(fc)exp(-Te,) 
q 

{B,{T)Ak) = -^*,(j)$,(fc)exp(-r6,) 



(2.31) 



whereas the equal-time contractions are given below 
.181). For the finite temperature contractions see c.f. 



F or lon gitudinal correlations the matrix-elements of S'f 
in (2.24) is given in a simple form for any position, I, 
therefore Gf jr) c an be obtained from the analogous ex- 
pression to ( |2.15| ) as 



G] 



q p>q 



■ exp[-T(eg + Cp) 



(2.32) 



III. PHENOMENOLOGICAL THEORY FROM 
SCALING OF RARE EVENTS 

In classical random ferromagnets where the critical be- 
havior is controlled by a random fixed point the distri- 
bution of several physical quantities (order-parameters, 
correlations, autocorrelations, etc.) is broad and as a 
consequence these quantities are not self-averaging: their 
average and most-probable or typical values are different. 
In random quantum spin chains the critical properties 
are expected to be controlled by the infinite-randomness 
critical fixed point ]3l| , where the distributions are ex- 
tremely (logarithmically) broad and as a consequence the 
average and typical behavior of these quantities are com- 
pletely different. The average is dominated by such re- 
alizations (the so called rare events), which have a very 
large contribution, but their fraction is vanishing in the 
thermodynamic limit. In this Section we identify these 
rare events for the random XY (and XX) model and 
use their properties to develop a phenomenological the- 
ory. Our basic observations are related to exact relations 
about the surface order-parameter and the energy of low- 
lying excitations. 



A. Surface order-parameter and the mapping to 
adsorbing random walks 

The local order-parameter at th e boundary is given by 
the simple formula in Eq. (2.21) as a sum of products 



of the ratio of the couplings J2j-i ^-nd Jf^. It is easy to 
see from Eq. (2.21) that in the thermodynamic limit the 



average surface order-parameter is zero (non-zero), if the 
geometrical mean of the Jf^ couplings is (smaller) grater 
than that of the ^Ij-i couplings. From this the definition 
of the control parameters in Eqs. ( p^ ) and (2.5) follows. 

Next we compute the average value of the surface 
order-parameter for the extreme binary distribution, i.e. 
the limit A ^ in ( ^.6[ ). For a random realiza- 
tion of the couplings the surface order-parameter at the 
critical point (ja = q = 1/2) is zero, whenever a prod- 
uct of the form of ni=i('^i^)~^j I — 1,2, . . . , L is infinite, 
i.e. the number of A-couplings exceeds the number of 
A^^-couplings in any of the [1,/] intervals. Otherwise 
the surface order-parameter has a finite value of 0{\). 
The distribution of the couplings can be represented 
by one-dimensional random walks that start at zero and 
make the i-th step upwards (for Jf^ = A~^) or downwards 
(for J|j = A). The ratio of walks representing a sample 
with finite surface order-parameter is given by the sur- 
vival probability of the walk Psurv, i-e. the probability of 
the walker to stay always above the starting point in L/2 
steps which is given by Psurv{L /2) ^ L^^l"^ . 

Next we consider the vicinity of the critical point, 
when the scaling behavior of the average surface order- 
parameter can be obtained from the survival probabilities 
of biased random walks [|l5| , where the probability that 
the walker makes a step towards the adsorbing bound- 
ary, q, is different from that of a step off the boundary, 
•p. The control parameter of the walk, 6^ = p — q, is 
analogous to the quantum control parameters Sa and 6d 
in Eqs. (2.2) and (p. 5]), respectively. Thus we have the 



basic correspondences between the average surface order- 
parameter of the XY (and XX) model and the surviving 
probability of adsorbing random walks: 

[mi{S,L)U ^ Psn.A5w,L/2), 6^5^, (3.1) 

We recall the asymptotic properties of the surviving 
probability of adsorbing random walks ||l5| . For unbi- 
ased walks: 

-Psurv(5» =0,L)~i-i/2^ (3.2) 
for walks with a drift away from the wall: 

Psui-viSw > 0, L ^ co) , (3.3) 
and for walks with a drift towards the wall: 



PsurviSw <0,L) ^ exp(-L/^u_,), 



in 



(3.4) 



In this way we have identified the rare events for the sur- 
face order-parameter, which are samples with a coupling 
distribution which have a surviving walk character. The 
scaling properties of the average surface order-parameter 
and the correlation length immediately follow from Eqs. 
( |3.2[ ), (3.3) and ( |3.4[ ) and will be evaluated in Section 
IV.A. 
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B. Scaling of low-energy excitations 



(3.8) 



The rare events controlling the surface order- 
parameter are also important for the low-energy excita- 
tions. Our results are obtained by using a simple relation 
for the smallest gap, ei(Z), of an open system of size /, 
i.e. with free boundary conditions, expecting that it goes 
to zero at least as ^ l/l. With this condition one can 
negle ct the r.h.s. of the eigenvalue problem of T in Eq. 
( 2.12 ), TVi — eiVi, and derive approximate expressions 
for the eigenfunctions $i and ^'i. With these one arrives 
at: 



mi m,] 



l/2-l 

-i-^f-i n 



jx 



(3.5) 



Here i 

parameter at the other end of the chain 



is defined in (2.21) and the surface order- 



IS given 

as in Eq. (|2^) replacing Jl^.i/Jfj by Jf+i_2]/ Ji-2j- 
(For details of the derivation of a similar expression for 
the quantum Ising chain see in Ref. P^.) 

Before using the relation in Eq. (|3.5|) we note that 
(surface) order and the presence of low-energy excitations 
are inherently related. These samples with an exponen- 
tially (in the system size) small gap have finite, 0(1), 
order-parameters at both boundaries and the coupling 
distribution follows a surviving walk picture. Such type 
of coupling configuration represents a strongly coupled 
domain (SCD), which at the critical point extends over 
the size of the system, L. In the off-critical situation, 
in the Griffiths phase the SCD-s have a smaller extent, 
I L, and they are localized both in the volume and 
near the surface of the system. The characteristic excita- 
tion energy of an SCD can be estimated from Eq. ( |3.5| ) 
as 



ei(0 



l/2-l 

n 



'J2j- 



'2j 



exp{-^ln(JVJ-)} , (3.6) 



where ?tr nieasures the size of transverse fluctuations of 
a surviving walk of length I and \n{J^ /Jv) is an aver- 
age ratio of the couplings, (it is ln(J'^/J°) for the XX 
model) . 

At the critical point {S = 0), where I ^ L, the size 
of transverse fluctuations of the couplings in the SC D i s 
Itr ~ ^1/2 |l|]. Consequently we obtain from Eq. ( |3.6| ) 
for the scaling relation of the gap: 



e((5 = 0, L) - exp(-const • L^/^) 



(3.7) 



Then the appropriate scaling variable is Ine/VX and the 
distribution of the excitation energy is extremely (loga- 
rithmically) broad. 

In the Griffiths phase the size of a SCD can be esti- 
mated along the lines of Ref. ||l^ as Z ~ In i and the 
size of transverse fluctuations is now /tr ^ ^ ~ In L. Set- 



ting this estimate into Eq. (B.6) we obtain for the scaling 
relation of the gap: 



where z is the dynamical exponent as defined in Eq. 



(1.1). The distribution of low-energy excitations can be 
obtained from the observation that an SCD can be local- 
ized at any site of the chain, thus Pl{1) ^ -PL(liie) ^ L. 
For a given large L the scaling combination from Eq. 
(|3.8D is Le^/^, thus we have: 



P(e) 



-l+l/z 



(3.9) 



As already mentioned z is a continuous function of the 
quantum control parameter 6 and we are going to calcu- 
late its exact value in Section V. 



C. Scaling theory of correlations 

The scaling behavior of critical average correlations 
is also inherently connected to the properties of rare 
events. Here the quantity of interest is the probabil- 
ity P^(Z), which measures the fraction of rare events of 
the local order-parameter mf. For the surface order- 
parameter mf it is given by the s urv iving probability, 



P^(l) = Psurv, according to Eq. (B.l). We start with 



the equal-time correlations in Eq. (2.3). In a given sam- 
ple there should be local order at both reference points of 
the correlation function in order to have C^(r) — 0(1). 
This is equivalent of having two SCD-s in the sample 
which occur with a probability of P2(l,l + r), which 
factorizes for large separation lim^^oo P2'(^ ^ + r) = 
P'^{l)P'^{L + r), since the disorder is uncorrelated. The 
probability of the occurrence of a SCD at position I, 
P^{1), has the same scaling behavior as the local order- 



parameter [mfjav = [{(t)p.\Si 
bulk point, < l/L < 1, as: 



K(i)]a 



which behaves at a 



(3.10) 



whereas for a boundary point, I ~ I, this relation in- 
volves the surface scaling dimension x^. Consequently 
Pf'il) transforms as P^(Z/6) = h''''^ P^{1) under a scal- 
ing transformation, when lengths are rescaled by a factor 
6 > 1 . Recalling that for spatial correlations there should 
be two independent SCD-s we obtain the transformation 
law: 



-22; 



(3.11) 



Now tak ing b = r one recovers the power-law decay in 
Eq. (2^) with the exponent 



(3.12) 



For critical time-dependent correlations the scaling be- 
havior is different from that in Eq. (3.11). This is due 
to the fact that disorder in the time direction is perfectly 
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correlated and the autocorrelation function in a given 
sample is Gi{t) = 0(1), if there is one SCD localised at 
position I. Therefore the average autocorrelation func- 
tion [G{'(lnr)]av scales as the probability of rare events 
P^'{l)■. 



[Gr(lnT)]e 



5--"[Gr(lnr/6i/2)]3 



(3.13) 



where we have used the relation in Eq. (4.3) between 
relevant length and time at the critical point. Taking the 
length scale as 6 = (Inr)^ we obtain for points / in the 
volume: 



[Gr(r)]av^(lnT)- 



(3.14) 



whereas for surface spins, I — 1, one should use the cor- 
responding surface decay exponent ry^. 

Next we turn to study the scaling properties of the 
average correlation functions in the Griffiths phase, i.e. 
outside the critical point. For equal-time correlations in 
a sample C^{r) = 0(1), if the SCD ext ends over a large 
distance of r, which according to Eq. (3.4) is exponen- 
tially improbable. Thus the average spatial correlations 
decay as 



[^^(r)],, ~ exp(-r/0, C ^ C« 



(3.15) 



where is defined in Eq. ( |3.4| ). On the other hand the 
autocorrelation function in a sample is G^(t) = 0(1), 
if there is one SCD localized at /, which occurs with a 
probability of P^{1) ^ 1/L. Consequently the average 
autocorrelation function, which scales as /'''(Oi trans- 
forms under a scaling transformation as: 

[Gr(r)]av = 6-MGr/b(T/6^)]av <5>0, (3.16) 

where we used the scaling combination t/V- in accor- 
dance with Eq. (1.1). Now taking h — r^/^ we obtain 



(3.17) 



for any type of autocorrelations, both in the volume and 
at the surface. 



SCD-s, having a coupling distribution of surviving RW 
character. The typical size of an SCD, as given by 
in Eq. (^^), is related to the average correlation length 
of the sys te m, [ ^jay. Then using the correspondences in 
Eqs. (|3.lD, (p]j) and (3.15) we get the relation: 



Klav - \S\- 



2 . 



(4.1) 



The typical correlation length, ^typ, as measured by 
the average of the logarithm of the correlation func- 
tion is different from the average correlation length. 
One can estima te th e typical value by studying the for- 
mula in Eq. 



( |2.21 



for the surface order-parameter, 
where the products are typically of Y\j{J 2j-i I J 2j)^ ^ 
exp(const • \5\L), thus [ms(L, (5(0)]typ ~ exp(— const • 
\5\L) ~ exp(— L/^typ)- Thus we obtain: 



Vtyp = 1 



(4.2) 



We note that at the critical point the largest value of 
the above products is typically of (•^Ij-i/^fj)^ ^ 
exp(^L^/^), since the transverse fluctuations in the cou- 
plings are of 0(L^/'^), thus we have [ms{L^6 = 0)]typ ~ 
exp(— const • L^/^). 



As shown in Eq. (3_^) the value of the smallest gap is 
related to the size of transverse fluctuations of an SCD, 
Ztr. Away from the critical point, when the correlation 
length is finite, one has ^tr ~ C^^^j and therefore the typ- 
ical relaxation time of a sample with typical correlation 
length ^ scales as 



In tr 



1/2 



(4.3) 



We note that the results in this part about length- 
and time-scales are valid both for the XY and XX mod- 
els. They also hold in identical form for the random 
TIM 1^,^, which can be understood as a consequence 
of the mapping of the ATy-chain into decoupled TIM-s. 
(See Appendix.) Since the corresponding scaling expres- 
sions for the random TIM have been studied in detail in 
previous numerical work ||r^ , p^ we do not repeat these 
calculations here. 



IV. CRITICAL PROPERTIES 



B. Quasi-long-range-order 



Here we consider in detail the random XY and XX 
chains in the vicinity of the critical points, as defined in 
Eqs. (2^) and (2^), respectively. The off-critical prop- 
erties of the systems in the Griffiths phase are presented 
afterwards in the following Section. 



A. Length- and time-scales 

As we argued in the previous Chapter the average be- 
havior of random quantum spin chains are inherently 
related to the properties of the rare events, which are 



At the critical point of random quantum chains the 
equal-time correlations decay with a power-law, (see Eq. 
( ^.4|) ), thus there is QLRO in the system. The decay ex- 
ponent of critical correlations are related to the scaling 
exponent of th e fra ction of rare events of the given 
quantity (see Eq. ( |3.12| )) and its value generally depends 
on the type of correlations of the disorder, thus it could 
be different for the XY and the XX models. Analyz- 
ing the scaling properties of the rare events in the XY 
and XX chains we have calculated the critical decay ex- 
ponents of different correlation functions, both between 
two spins in the volume and for end-to-end correlations. 
Our results are presented in Table I. 
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ri'-{XY) 


vHxx) 


bulk 




2(*) 


4 


2(*) 


surface 


1 


1 


2 


1 



TABLE I: Decay exponents of critical correlations in the ran- 
dom XY and XX chains. The exponents with a superscript 
are those calculated by Fisher with the RG method y, 
whereas '**■' follows from the results of the random TIM in 
Ref. §. 

In the following we are going to derive these expo- 
nents by analytical and scaling methods and then con- 
front them with the results of numerical calculations. 

1. Longitudinal order-parameter 



as shown in Table I. 

We studied the order-parameter profile [?7if]av numeri- 
cally for large finite systems up to L = 256. As shown in 
Fig. H the numerical points of the scaled variable L[m^]av 
are on one scaling curve rh^{l) for different values of L. 
The scaling curve has two symmetric branches for odd 
and even lattice sites, which cross at I = L/2. The up- 
per part of the curves in the large L limit is very well 
described by the function rh^{l)^ = ^sin(7rZ)^^/^, which 
corresponds to the conformal result about off-diagonal 
matrix-element profiles p^ ]: 

HU-{jj (sin^-j , (4.6) 



We start with the scaling behavior of the longitudinal 
order-parameter m^, which in t he XX chain is given by 
the simple formula in Eq. ( ^.26| ) . Summing over all sites 
one gets the sum-rule 



1=1 



= 1 XX - model , 



(4.4) 



where we have used Eq. (2.13) and the fact that the <i>^ 
and are normalized. Since this sum-rule is valid for 
the average quantities, too, we get immediately 



H Jav 



= L- 



"(l/L) , 



(4.5) 



where m^(Z) is a scaling function with I = l/L. Conse- 
quently for bulk spins the finite-size dependence of the 
l ocal order-parameter is [fTi^Jav ^ L~^, thus from Eq. 
( |3.10D we have x%XX) = 1 and from Eq. (|3.12|) the 
decay exponent is 

r]'{XX) = 2 

as given in Table I. A further consequence of the sum-rule 
in Eq. (4.4) is that the average value of the bulk order- 
parameter is the same, if the averaging is performed over 
any single sample. Thus the order-parameter and the 
correlation function {S[S[^^) are self-averaging. This is 
quite special in disordered systems where the correlations 
are generally not self-averaging [Q. The self-averaging 
properties of the S'^-correlations provides an explanation 
of the accurate numerical determination of the decay ex- 
ponent r]^{XX) = 2 in previous numerical work [ p5|j2^ . 

The surface order-parameter mf for the XX model sat- 
isfies the relation mf = 2 (m^)^, which follows from Eqs. 
( ^.2l[) and (|2.26|) . Then a rare event with mf = 0(1) 
is also a rare event for the order-parameter mf , conse- 
quently the fraction of rare- even ts is given by the 
surviving probability in Eq. (3.2). Thus the scaling di- 
mension is xf = 1/2 and the decay exponent of critical 
end-to-end correlations is 

vKxx) = 1 



with = 1 and xf = 1/2. On the other hand the lower 
part of the curves in Fig. ^ is given by m^(i)\ = ^sin(7rf), 
which corresponds to Eq. ( |4.6| ) with a;| = 2. Thus 
we obtain that average critical correlations between two 
spins which are next to the surface are decaying as 
[C^(2, L — l )]av ~ Using the sum-rule for the pro- 

file in Eq. (4.4) and the conformal predictions one can 
determine the pre-factor A from normalization. Then 
from the equation A/2 ^^[{smTix)~^^'^ -\- sinTrxJda; = 1, 
one gets A = .86735, which fits well the numerical data 
on Fig. |. 
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ra 6 
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0.867 rsin(7rx)]"^'^ 
0.867 sin(7rx) 
L=8 
L=16 
L=32 
L=64 
L=128 




(l-0.5)/L 

FIG. 2. Finite size scaling plot of the longitudinal or- 
der-parameter profiles [mf]av for the XX model at critical- 
ity for different system siz es ca lculated numerically with the 
fermion method using eq. {2.2t). The data are for the uniform 
distribution, averaged over 50000 samples. The conformal re- 
sults are indicated by full lines. 

These results about the conformal properties of the 
profile are in agreement with similar studies of the ran- 
dom TIM [|l^,^. Thus it seems to be a general trend 
that critical order-parameter profiles of random quantum 
spin chains are described by the results of conformal in- 
variance, alth ough these systems are strongly anisotropic 
(see Eq. (4.3)) and therefore not conformally invariant. 
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Next we turn to study the order-parameter mf and 
the longitudinal correlation function in the random XY 
model. In this model the disorder in the Jf and Jf cou- 
plings is uncorrelated, therefore one can perform aver- 
aging in the two subspaces To- and Tr, or in the two 
decoupled TIM-s, independently. Note that the expres- 



sion for mf in Eq. ( 2.25| ) is given as a product of two 



vector-components, where each vector belongs to differ- 
ent subspaces and have the same average behavior. Since 
the couplings entering the two separate eigenvalue prob- 
lems are independent one gets for the disorder average 



[mf]av = [*l(/)]av • [*2(0]a 



(4.7) 



Since the probability for mf being of order one is the 
product of the probabilities for and 5*2(0 being of 

order one we conclude that the scaling dimension for mf 
in the random XY chain is twice that for the random 
XX chain. Thus the decay exponents are 



and 



VKXY)^2 



in the bulk and at the surface, respectively, as shown in 
Table 1. 

The numerical results about the order parameter pro- 
file is shown in Fig. ||. The data collapse is satisfactory, 
although not as good as for the XX model. Similar con- 
clusion holds for the relation with the conformally pre- 
dicted profile, which is also presented in Fig. |^. 




(l-0.5)/L 

FIG. 3. Finite size scaling plot of the longitudinal or- 
der-parameter profiles [mf]av for the XY model at critical- 
ity for different system size s calc ulated numerically with the 
fermion method using eq. ( 2.26 ). The data are for the uni- 
form distribution, averaged over 50000 samples. 



its average behavior follows from the adsorbing random 
walk mapping in Section UFA. Then from Eqs. ( |3.l|) and 
.21) one gets xf = 1/2 and 



1 



both for the random XY and XX models, as shown in 
Table F The value of the decay exponents follows also 
from the mapping to two TIM-s. As shown in Eq. ( [A^ ) 
in the Appendix the correlation function {S21S', 



IS 



expressed as the product of spin correlations in the two 
TIM's, one with open boundary conditions, but the other 
is taken with fixed-spin boundary conditions in terms of 
dual variables. For end-to-end correlations this second 
factor in the product is unity, since it is the correlation 
between two fixed spins. Therefore end-to-end correla- 
tions between the random TIM and the random XY and 
XX models are identical and the decay exponent corre- 
sponds to the value in Table I. 



0.45 



CO 

o 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
(I-0.5)/L 

FIG. 4. Transverse order-parameter profile [mf]av for the 
XY model at criticality for different system size s calc ulated 
numerically with the fermion method using eq. (2.22). The 



data are for the uniform distribution, averaged over 50000 
samples. 

For bulk correlations one can easily find the a nswe r 
for the XY model with the mapping in Eq. (|A^). 
When the two points of reference are located far from 
the boundary the boundary condition does not matter 
and after performing the independent averaging for the 
two factors of the product one obtains [('S'2;5'2j_|_2r)]av = 



l/4[ 



I "i+r/Jav 



thus 



'n'^iXY) = 2r]{TIM) = 3 - Vs , 



(4.8) 



2. Transverse order-parameter 



We start with the surface order-parameter. 



as 



given by the simple formula in Eq. (2.21). This for- 
mula is identical both for the XY and XX models and 



where the last result follows from Fisher's RG calcula- 
tion 1^. (As shown in Ref. js^] the rare events for the 
bulk order-parameter in the TIM are samples having a 
coupling distribution of average persistence character). 
The scaling exponent (XY) can identically be obtained 
from the expression of the order-parameter profile in Eq. 
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(|^, which is in the form of a product of the two Ising .|__^-| 
order-parameters and for the XY model the two factors 
are averaged independently. 

For the XY model the numerically calculated profile is 
shown in Fig. ^. The scaling plot with the exponents in 
Table I is reasonable, although larger systems and even 
more samples would be needed to reach the expected 
asymptotic behavior, as predicted by conformal invari- 
ance in Eq. (4.6). 

The arguments leading to the prediction ( |4.8| ) for the 
transverse bulk order parameter exponent do not apply 
for the XX model and one cannot obtain a si mple es- 
timate for the bulk decay exponent from Eqs. (A8) or 
(2.22) due to the following reason. The expressions with 
the parameters of the two quantum Ising chains contain 
real and dual variables for the two (cr and t) systems. 
Since Jf = Jf — Ji a. domain of strong couphngs in the 
a chain corresponds to a domain of weak couplings in 
the T chain and reverse. Therefore the rare events of the 
TIM can not be simply related to the rare-events of the 
XX chain. 

The value for r]^{XX), however, can be obtained by 
the following argument. For simplicity let us consider 
the extreme binary distribution in which J2i = 1 and 
J2i_i = A or 1/A with pro babil ity 1/2, taking the limit 
A ^ 0. Then, from Eq.( 2.21 ), one gets only then a 
non-vanishing transversal surface magnetization, when 
the disorder configuration has a surviving walk character 
(meaning Hi^i J^i-i < oo for alH = 1, . . . , L/2 — 1). 
This implies, also for general distributions of couplings 
that mf ~ 0(1) only if the surface spin is weakly coupled 
to the rest of the system. It is instructive to note the dif- 
ference to the surface magnetization in the TIM, where 
mf ~ 0(1) when the surface spin is strongly coupled to 

the rest of the system, meaning that 11^=1 (iZ-^i) < °^ f^^' 
all Z = 1, . . . , L — 1 for the extreme binary distribution. 

The same remains true for a bulk spin, which also 
has non-vanishing transverse magnetization only if it is 
weakly coupled to the rest of the system (the trivial ex- 
ample being when both its couplings to the left and to 
the right are exactly zero, which gives the maximum 
value mf — 1/2). Thus the central spin in a chain of 
length, say 2L — 1, has ~ C>{1) if and only if the 
bond-configurations on both sides have surviving char- 
acter, as it is depicted in fig. ^ for the extreme binary L"'l/2J 
distribution. Since the probability Psurv{L/2) for a con- 
figuration of L /2 couplings to represent a surviving walk 
is Psurv{L/2) - L-^/"^ it is 

mf ^ {P,urv{L/2)Y - L~\ i.e. x^XX) = 1 . (4.9) 



-5 -3 -1 q 1 3 5 



L-1 



From this one obtains 

?7^(XX) = 2 , (4.10) 
as given in Table I. 



Ill Vk X X 111 X l/X X l/X 

1 1 \ \ ^ 1 \ 1 1 



FIG. 5. Sketch of a bond configuration for a chain of length 
2L — 1 that gives a non-vanishing transverse magnetization 
0{1) for the central (bulk) spin. The example is for 
the extreme binary distribution. Weak couplings (J2i-i = A) 
correspond to downward steps of the random walk on both 
sides of the central spin (here at 0). Note that both, the right 
and the left half of the random walk have surviving character, 
i.e. do not cross the starting point. 

We verified the strong correlation between weak cou- 
pling and non-vanishing transverse order parameter nu- 
merically in the following way: We considered a chain 
with L + 1 sites and the couplings at both sides of the 
central spin were taken randomly from a distribution 
called SW iQ, which represents those samples in the 
uniform distribution which has a surface magnetization 
of mf{SW) > 1/4. (Thus cutting one of the couplings 
to the central spin results a local magnetization greater 
than 0.25.) Then we calculated numerically the order- 
parameter at the central spin and its average value over 
the SW configurations [m^/jlsw as given in Table I. 



L 


2[mf]sw 




16 


0.817 


0.531 


32 


0.806 


0.471 


64 


0.799 


0.431 


128 


0.792 


0.413 


256 


0.791 


0.383 



TABLE II: Surface and bulk transverse order-parameters av- 
eraged over 50000 SW configurations for the uniform distri- 
bution. 

As seen in the Table the averaged surface order- 
parameter stays constant for large values of L, whereas 
the bulk order-parameter decreases very slowly, actu- 
ally slower than any power. The data can be fitted by 
(InL)""', with a w 0.5. Thus we conclude 
t hat t he numerical results confirm the exponents given in 
( 4.10 ), however there are strong logarithmic corrections, 
which imply for the average transverse correlations 



ln"^(r) XX -model 



(4.11) 



These strong logarithmic corrections render the numeri- 
cal calculation of critical exponents very difficult |^,^ . 
In earlier numerical work using smaller finite systems dis- 
order dependent exponents were reported | p5[ |. We be- 
lieve that these numerical results can be interpreted as 
effective, size-dependent exponents and the asym ptotic 
critical behavior is indeed described by Eq. ( 4.11 ) . 
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Note that our results in Table I satisfy the relation 
r]^{XX) = ri^{XX), both in the volume and at the sur- 
face, which corresponds to Fisher's RG result In 
this way we have presented in independent justification 
of Fisher's RS phase picture, where the average corre- 
lations are dominated by random singlets, so that the 
distance between the pairs could be arbitrarily large. 



is clear that x^{L) grows at least up to the theoretical 
limit cc^ = 1, although it could, in principle, reach even a 
larger value. We note that similar observation was made 
by Henelius and Girvin from the average S'^ correlation 
function, where the effective rj^ exponents seem to grow 
over the theoretically predicted value oi rj^ — 2. (See 
Fig. 2 of Ref 




C. Autocorrelations 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



(l-0.5)/L 

FIG. 6. Transverse order-parameter profile [mf]av for the 
XX model at criticality for different system size s calc ulated 
numerically with the fermion method using eq. (2.21). The 
data are for the uniform distribution, averaged over 50000 
samples. 



We checked numerically the above theoretical predic- 
tions in the random XX model. In Fig. ^ we present 
the scaled mf profiles for the binary distribution for fi- 
nite systems up to L = 512. The profiles have a broad 
plateau and the points of mf do not perfectly fall on 
one scaling curve due to strong finite-size effects. Even 
system sizes as large as L = 512 appear to be insufficient 
to get rid of such correction terms. Therefore we have 
calculated the effective size-dependent exponents 
by a two-point fit. For this we have averaged the order- 
parameter in the middle of the profile for L/i < I < 3L/A 
and compared this average values for finite systems with 
L/2 and L sites. As seen in Table III the effective ex- 
ponents are monotonously increasing with the size of 
the system and they are not going to saturate, even for 
L = 512 M. 



L 




16 


0.635 


32 


0.677 


64 


0.730 


128 


0.823 


256 


0.872 


512 


0.910 



TABLE III: Effective bulk scaling dimension of the transverse 
order-parameter in the random XX chain. 

From the data in Table III one can not make an ac- 
curate estimate about the limiting value of (L) , but it 




6 8 
In T 

FIG. 7. Spin autocorrelation function [Gj'(r)]av for the 
XX model for L = 32, 64 and 128 calcu lated nume rically 
with the fermion method using eqs. ( 2.27 ) and ( 2.32 ). The 
data are for the binary distribution (A — 4), averaged over 
50000 samples, a) (Top) shows 1 = 1, the surface transverse 
autocorrelations, b) (Bottom) shows I = L/2, the bulk longi- 
tudinal autocorrelations. 

According to the scaling theory in Section III.C the de- 
cay of average critical autocorrelations in random quan- 
tum spin chains is ultra-slow, it take s place in logarith- 
mic time-scales, as given in Eq. (3.14). Here we confront 
these predictions with the results of numerical calcula- 
tions. We start with the surface autocorrelation function 
[Gf{T)]av for the XX model, which is calculated in the bi- 
nary distribution (A = 4) on finite systems up to L = 128. 
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As seen in fig. |^ (top) the logarithmic time-dependence is 
well satisfied and the decay exponent is found in agree- 
ment w ith rj T(XX) = 1 as given by the scaling result 
in Eq. ( 3.14 ). For bulk spin critical autocorrelations we 
considered [G^^2(''')]av for the XX model. Again the 
numerical results in Fig. |^ (bottom) are consistent with 
a logarithmic decay with an exponent ri^-{XX) = 2, as 
given in Table I. 




x= 10 
T= 100 
T= 1000 
T= 10000 
T= 100000 



^ 0.1 



0.01 



1 2 3 4 5 6 

X = -In Gi(t) / Inx 



T= 10 
T= 100 

T= 1000 
T= 10000 
T= 100000 



Q- 
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X = -In Gf/Jx) / Inx 
FIG. 8. Scaling plot of the probability distribution of the 
autocorrelation function Gf(r) for the XX-model for differ- 
ent values of r at criticality (L = 128). The data are for the 
uniform distribution averaged over 100000 samples, a) (Top) 
shows / = 1, the surface transverse autocorrelations, b) (Bot- 
tom) shows I = L/2, the bulk longitudinal autocorrelations. 



Next we turn to study the distribution of critical au- 
tocorrelations. As we have seen the average behavior 
is logarithmically slow, but for typical samples, as de- 
scribed in Appendix B, one expects a faster decay with 
a power-law time-dependence. Then G'['(r) ~ t""*" and 
the 7 exponent could vary from sample to sample. Such 
type of "multi-scaling" behavior of the autocorrelations 
has been recently observed by Kisker and Young |38| in 
the random quantum Ising model. In Fig. |^ we have nu- 
merically checked this assumption for the critical auto- 
correlations G'i(t) and G|^2(''')j respectively, of the ran- 



dom XX chain, the average behavior of those have been 
studied before. As seen in Fig. || we have obtained in- 
deed a good data collapse of the probability distributions 
P''(7) in terms of the scaling variable 7 — — InCf'/lnr 
for both type of autocorrelations, but the scaling curve 
in the two cases are different. 

The average correlation function generally have a con- 
tributions from the scaling function, P^(7), but there 
could be also non-scaling contributions, as found for the 
random quantum Ising chain in Ref p9| . The scaling 
contribution is coming from the small 7 part of the scal- 
ing function, which according to Fig. || (top) for the au- 
tocorrelations G'i(t) approaches a finite value linearly, 
P^il) ^ A + B-f. Thus we have for the average autocor- 
relations: 



[GUr)U = / P^(7)GT(T)d7 
Jo 



(A + B"f) cxp(— 7 In r)(i7 



A(ln^ 



(4.12) 



in agreement with the scaling result in Eq. (3.14) and 
with the numerical result in Fig. || (top) We note that 
the correction to scali ng co ntribution to the average au- 
tocorrelations in Eq. ( 4.12| ) is also logarithmic. 



For the critical autocorrelation G|^^2(''') the scaling 
function in Fig. |^ (bottom) for small 7 approaches lin- 
early zero [|oj P^i'y) ~ 7. Thus the scaling contribution 
to the avera ge au tocorrelation, as evaluated along the 



lines of Eq. ( 4.12 ), is [G'f(T)] ay ^ (I nr) ^, in agreement 



with the scaling result in Eq. (3.14) 



V. GRIFFITHS PHASE 

Random quantum systems exhibit unusual off-critical 
properties: they are gap- less in a extended region, < 
\S\ < Sg, as a result of the so called Griffiths-McCoy 
singularities ||l^,^. In this Griffiths phase the system 
is critical in the time direction, although spatial correla- 
tions decay exponentially. 

Quantitatively the basic information is contained in 
the dist ribu tion of low energy excitations, P{e), as given 
in Eq. (3.9). With this the average autocorrelations can 
be obtained as: 



P{e) exp(— Te)(ie ~ t 



-l/z 



(5.1) 



which is expected to hold for any component of the spin 
p2[ . In th is way we have recovered the scaling result in 
Eq. ( 3.17 ). In the Griffiths phase also some thermody- 
namic quantities arc singular, which are expressed as an 
integral of the autocorrelation function. We mention the 
local susceptibility xf at site /, which is de fined through 
the local order-parameter mf in Eq. ( 2.16| ) as 

^f= lim (5.2) 
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where Hf is the strength of the l ocal longitudinal field, 
which enters the Hamiltonian in {2A_) via an additional 
term HfSf. xf can be expressed as: 



Xt 



2E 

("I 



i(»i^rio)i 

E„ — En 



(5.3) 



thus its average value scales in finite systems as xf{L) ^ 
L^~^, where we have used the scaling relation in Eg . 
and the fact that the matrix-element in Eq. ( |5.3| ) 
1/L, since an SCD can be localized at any site of 
For a small finite temperature T we can use 
e and we have for the 



(3 



IS 

the chain 

the scaling relations T 
singular behavior: 



T 



-l+l/z 



(5.4) 



To estimate the temperature dependence of the aver- 
age specific heat, [C(T)]av, we calculate first th e av er- 
age excitation energy per SCD with P{e) in Eq. ( |3.9| ) as 
J eP{e)de ^ e^/^"*"^, which is proportional to the thermal 

rpi/z+i fj-Qjj^ which we obtain: 



excess energy per spin 
[C(T) 



(5.5) 



We note that several other physical quantities are singu- 
lar in the Griffiths phase (non- linear susceptibility, higher 
excitations, etc) and the corresponding singularities are 
expected to be related to the dynamical exponent z. For 
a detailed study of this subject in the random quantum 
Ising model see Ref. 
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In the following we calculate the exact value of the dy- 
namical exponent using the same strategy as for the ran- 
dom quantum Ising model in Ref |po| , ^ . Our basic ob- 
servation is the fact that the eigenvalue problem of the To- 
(or Tr) matrix can be mapped through an uniter transfor- 
mation to a Fokker-Planck operator, which appear in the 
Master equation of a Sinai diffusion, i.e. random walk in 
a random environment p^ . The transition probabilities 
of the latter problem are then expressed with the cou- 
pling constants of the spin model. The Griffiths phase 
of the spin model corresponds to the anomalous diffusion 
region of the Sinai walk and from the exact results about 
the scaling form of the energy scales in this problem one 
obtains for the dynamical exponent of the XY model: 



Jy 



l/z- 



(5.6) 



whereas for the XX model the result follows with the 



corresp onde nces in Eq.(2.8). For the binary distribution 
in Eq. (2.6) the Griffiths phase is for 1 < Jq < A and z 
is given by: 



cosh 



In A 

z 



(5.7) 



Tn(l-z-2) =-lnJo^ , 



(5.8) 



and the Griffiths phase extends to 1 < Jq < cx). 

Next we are going to study numerically the Griffiths 
phase and to verify some of the scaling results described 
above. In this respect we shall not consider those quanti- 
ties which have an equivalent counterpart in the random 
quantum Ising model (distribution of energy gaps, lo- 
cal susceptibility, specific heat, etc), since that model 
has already been thoroughly investigated numerically 
The autocorrelation functions, however. 



|l3|,^|l5y2^]. 

are different in the two models and we are going to study 
those in the following. 
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FIG. 9. The average surface (top) and bulk (bottom) au- 
tocorrelation function [G^y2 i(''')]av of the XX model in the 
Griffiths-phase for various values of ho. The straight lines 
have a slope of l/z{ho), where the dynamical exponent zjho ) 
agrees well the exact value determined via the formula (5.6). 



For the uniform distribution 



The data are for the uniform distribution averaged over 50000 
samples of size L — 128. 

The average bulk longitudinal autocorrelation function 
[^L/2(''')]av of the XX model is shown in Fig. || in a log- 
log plot at different points of the Griffiths phase. The 
asymptotic behavior in Eq. (5.1) is well satisfied and 
the dynamical exponents obtained from the slope of the 
curves are in good agreement with the analytical results 
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in Eq. (5.6). Similar conclusion can be drawn from the 
average surface transverse autocorrelations, [Gf(r)]av, as 
shown in Fig. ||. 
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FIG. 10. (Top): Probability distribution of the bulk lon- 
gitudinal autocorrelation function of the XX model in the 
Griffiths phase for ho = 1.5. The data are for the uniform 
distribution from 100000 samples of size L — 128. 
(Bottom): Scaling plot of the data in the top figure. The 
scaling variable [In G(t)]/t^^'^'''^^ contai ns t he dynamical ex- 
ponent z{ho) known from the formula (5.6). The full curve 
is the theoretical prediction in ( 5.11 ) using the exact value of 
z{ho — 1.5) — 2.659 and a fit-parameter c = 0.22. 

Next we study the distribution of the autocorrelation 
functions. In Fig. 10 the distribution of the bulk longitu- 
dinal autocorrelation function of the XX model is given 
at different times t. As argued in the Appendix the typ- 
ical autocorrelations are of a stretched exponential form 



G(r) - exp (-const t^/(^+^^ 

thus the relevant scaling variable is 
In G(t) 



-1/(2 + 1) 



(5.9) 



(5.10) 



Using this scaling argument we obtained a good data col- 
lapse of the points of the distribution function as shown 



in Fig. . We note that for the random quantum Ising 
model Young |l^ has also derived the scaling function 
from phenomenological arguments, 



P{x) = c(ca;)^/^ exp 



1 



-{cx) 



(5.11) 



which is also presented in Fig. ^ One can see consid- 
erable differences between the numerical and theoretical 
curves. Similar tendencies have been noticed for the ran- 
dom quantum Ising model in Ref jl6j . The discrepancies 
are probably due to strong correction to scaling or finite 
size effects. These corrections, however does not affect 
the scaling form in Eq. 



(5.10) 



VI. DISCUSSION 

In this section we first discuss the possible extension of 
our results to random XX Z chains and to higher dimen- 
sional systems, then we conclude with a brief summary 
of our findings. 



A. Random XXZ chains 



The more general XXZ (or XYZ) Heisenberg spin 
chain, where the Hamiltonian in Eq. (2T) contains an 
additional interaction term of the form 



L-l 



QZ QZ 



(6.1) 



1=1 



can be treated perturbatively when <C J^,J^. Here 
we consider the XXZ chain, with Jf — XJi and A ^ 1. 
To see the scaling behavior of the energy gap we express 
the small perturbation, Hz, in terms of two decoupled 
TIM-s (see Appendix) as 

L/2-1 

^^ = -i E (•^2Vi<T-f -f J>f af+iTf rf+O , (6.2) 
1=1 

whose expectation value in the unperturbed ground state 
is just the product of the local energy-densities of the 
TIM-s. The perturbative correction to the gap, Az{L) = 
{l\Hz\l) — (Olff^lO), evaluated with the states of the un- 
perturbed Hamiltonian, (0| and (1|, is proportional to the 
gap of one of the random TIM's, ea{L). At the critical 
point ^aiL) ^ exp(— const • L^^^) |jl^, which is the same 
scaling form as for the X Y m odel in Eq.(3.7). Thus the 
scaling relation in Eq. (O) is valid also for the XXZ 
chain. In the Griffiths-phase one has again e(j(i) ^ L~^, 
and the dynamical exponent z, has the same value as in 
Eq. (5.6). Thus we arrive at the conclusion that also in 



the Griffiths phase the corresponding scaling relation in 
Eq. (1.1) is valid in the same form for the random XXZ 
chain, at least for small couplings. 
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Next we study the asymptotic properties of average 
critical correlations in the random XX Z model through 
the scaling behavior of the local order-parameters. As 
we argued in Section III these quantities are related to 
the fraction of rare events, P'^, and here we are going to 
investigate the influence of the perturbation Hz to P^. 
We start with the surface transverse order-parameter 
and recall that it is maximal, i.e. m^f = 1/2, if the surface 
spin is disconnected in the XF-plane, i.e. Jf = Jf — 0. 
Evidently the value of mf = 1/2 does not change for 
any finite value of the coupling Jf . Now consider the 
infinite-randomness fixed point of the XX chain with 
the extreme binary distribution, where a rare event is 
represented by couplings with a surviving random walk 
configuration and with m\ = 0(1). Roughly speaking, a 
rare event is formally equivalent to a situation, in which 
there is a very weak surface coupling of Ji = 0(L^^/^), 
where L is the system size. Then switching on homo- 
geneous and finite couplings the lowest excitation of 
the chain stays localized at the surface, since the shape of 
the wave-function does not change significantly in first or- 
der perturbation theory. Consequently the surface order- 
parameter is still rrii = 0(1), and the sample is a rare 
event for the XXZ chain, too. For small random cou- 
plings Jf the accumulated fluctuations in Jf are diver- 
gent as ~ L^/^, however these are still negligible com- 
pared with the fluctuations in the transverse couplings. 
Thus the rare events of the XX chain are identical with 
those appearing in the XXZ chain for small values of the 
random longitudinal couplings. As a consequence the 
critical end-to-end average correlations decay with the 
same exponent as given in Table I. Since the rare events 
for other local order-parameters are also connected to 
SCD-s with localized wave-functions the stability of the 
infinite-randomness fixed point holds for the other crit- 
ical correlations, too. Actually, it seems to be plausible 
that the attracting region of the XX fixed point extends 
up to [In J^]av > [In J^]av, i.e. where the average trans- 
verse couplings are larger than the longitudinal ones, thus 
up to the random XXX fixed point, in agreement with 
Fisher's conjecture [Q. 



B. Higher dimensions 

In one dimension the topology is special since there is 
only a single path between two points, whereas in higher 
dimensional lattices one has several distinct paths con- 
necting two points. This topological difference is essen- 
tial when random XX magnets are considered in higher 
dimensions. Let us consider again the surface transverse 
order-parameter and construct a rare event in the ex- 
treme binary distribution. For this purpose the surface 
spin should be extremely weakly coupled to the bulk of 
the system. Thus considering any non-self-crossing path 
from the spin to the volume one should have a surviving 
random walk configuration in the couplings. In higher di- 



mensions the number of such paths grows exponentially 
with the size of the system, L, thus the fraction of rare 
events, which is related to the length as a power in one di- 
mensions, becomes exponentially small in higher dimen- 
sions. Consequently the infinite-randomness fixed point 
picture is not applicable here and one concludes that the 
critical properties of higher dimensional XX and Heisen- 
berg antiferromagnets are controlled by conventional ran- 
dom fixed points. This result is also in agreement with 
numerical RG calculations in 2d []l2|. We note that in 
contrast to random Heisenberg antiferromagnets the ran- 
dom ferromagnetic quantum Ising models in higher di- 
mensions are still controlled by infinite-randomness fixed 
points [||,0. 



C. Summary 

Quantum spin chains in the presence of quenched dis- 
order show unusual critical properties, which are con- 
trolled by the infinite-randomness fixed point. A common 
feature of these systems is that various physical proper- 
ties, especially those related to local order-parameters 
and correlation functions are not self-averaging and their 
average behavior is determined by the rare events (or rare 
regions), which give the dominant contribution, although 
their fraction is vanishing in the thermodynamic limit. In 
this paper we have performed a detailed study of the scal- 
ing behavior of rare events appearing in the random XY 
and XX chains. We identified the rare events as strongly 
coupled domains, where the coupling distribution follows 
some surviving random walk character. From the scal- 
ing properties of the rare events we have identified the 
complete set of critical decay exponents and found exact 
results about the correlation length exponent and the 
scaling anisotropy. 

Another new aspect of our work was the study of dy- 
namical correlations. We have obtained the asymptotic 
behavior of the average autocorrelation function and de- 
termined the scaling form of the distribution of autocor- 
relations. In the off-critical regime we investigated the 
singular physical quantities in the Griffiths phase. In 
particular we have obtained exact expression for the dy- 
namical exponent z, which is a continuous function of 
the quantum control-parameter and the singularities of 
all physical quantities can be related to its value. 
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APPENDIX: MAPPING TO DECOUPLED ISING 
QUANTUM CHAINS 

We start here with the observation in Section IIB that 



the eigenvalue matrix T in Eq. (2.12) can be represented 
as a direct product T = ^T^.. The trigonal matri- 
ces Tg. , Tt of size L X L represent transfer matrices of 
directed walks, which are in one-to-one correspondence 
with Ising chains in transverse field |^ defined by the 
Hamiltonians: 



^ L/2-1 



1=1 

L/2-1 
i=l 



L/2 

4 ^ 

i 

L/2 

4 



^ LI2. 
1 - -^^^21-1' 



(Al) 



Here the trf and rf are two sets of Pauli matrices at 
site i and there are free boundary conditions for both 
chains. We can then write Hxy — Ha + -f^r- Note the 



symmetry cr. 



and Jf ^ Jf, thus anisotropy 



in the XY model has different effects in the two Ising 
chains. 

One can easily find the transformational relations be- 
tween the XY and Ising variables: 



2i-l 



2i-l 



whereas the inverse relations are the following: 



(A2) 



25: 



2i-l 



-'2i-l 



i-l 



z n Qx : 



II' 



n 

i 



(A3) 



We note that a relation between the XY model and two 
decoupled Ising quantum chains in the thermodynamic 
limit is known for some time here we have ex- 

tended this relation for finite chains with the appropri- 
ate boundary conditions. These are essential to map local 
order-parameters and end-to-end correlation functions. 
End-to-end correlations are related as 



L/2 



(A4) 



since in the ground state (Hiii '^i) = 1- Similarly 



thus the end-to-end correlations in the two models are in 
identical form. As a consequence the corresponding de- 
cay exponent in the random models, 77^ in Table I is the 
same in the two systems and the same conclusion hold s 
also for the correlation length exponent, v in Eq. (4.1). 
These results are also independent of the type of correla- 
tion of the disorder, thus are valid both for the XY and 
XX models. 

Correlations between two spins at general positions 21 
and 21 + 2r are related as 



1 



>l "l+r 



i+ii 



(A6) 



The second factor in the r.h.s., (ni=i '''i+i)' defines a 
string-like order-parameter |p6|| what can be expressed in 
a simpler form in terms of the dual Ising variables '?^]^/2' 
which are defined on the bonds of the original Ising chain 
as 



~Z n-X n-X 

^i+1/2 — M U 



+ 1 



Z 2LX ~X 

~ ^1-1/2^1+1/2 



(A7) 



Under the duality transformation fields and couplings are 
exchanged, therefore the vanishing bonds at the two ends 
of an open chain are transformed to vanishing fields, thus 
the dual chain has two end spins fixed to the same state. 
So we obtain for the correlations in Eq. (A6) 



/ QX QX \ 

W2;'-'2;+2r/ 



1 



cr, a, 



I ";+r)(^/+l/2'i+r+l/2/ 



(A8) 



where the superscript denotes fixed-spin boundary 
condition. For non-surface points the average value of 
the correlation function in Eq. (A8) depends on the type 
of disorder correlations. For the XY model, where the 



disorder is uncorrelated the two factors in Eq. (A8) can 
be averaged separately, whereas this is not possible for 
the XX model. We treated this point in Sections IV. B. 2. 



APPENDIX: DISTRIBUTION OF 
AUTOCORRELATION FUNCTIONS 



The autocorrelation functions is represented by the 
general form: 



G(t) 



E 



A4pexp(-TA£:fe) 



(Bl) 



where the dominant contributions to the sum in Eq. (Bl) 
are from SCD-s, which are localized at some distance I 
from the spin and have a very small excitation energy, 
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AE{1). The scaling form of AE{1) follows from the con- 
siderations in Section III.B and one obtains from Eqs. 
(EtI) and (K~ 



AE{1) 



(B2) 



at the critical point and in the Griffiths phase, respec- 
tively, where eo denotes the energy scale. Thus the larger 
the distance from the spin the larger the probability to 
have an SCD with a very small energy. For the matrix- 
element, |M(/)p, the tendency is the opposite since the 
overlap with the wave-function of the SCD is (exponen- 
tially) decreasing with the distance. The corresponding 
scaling form can be read from the typical behavior of the 
surface order-parameter as given below and above Eq. 
(U) as 



|M(0|^ 



exp(- 
exp(- 



,5 = 

-V6yp), S<0 



(B3) 



Then G{t) in Eq. (Bl) can be approximated by a sum 
which runs over SCD-s localized at different distances, 
I, and this sum is dominated by the largest term with 
l = lo- 



G(r) ~ |M(;o)|'exp(-rAS(/o)) 



(B4) 



Using the scaling forms in Eqs. ( |B2[ ) and (B3) one gets 
following result. 

At the critical point the characteristic distance is = 
[\n{TeoA/B)/A]'^ and the typical autocorrelation function 
decays as a power: 

G(T)-r-^/^, 6^0, (B5) 

Thus the relevant scaling variable of the problem is 

In G(r) 



7 = 



In^ 



5 = 



(B6) 



In the Griffiths phase the characteristic distance has a 
power-law t dependence, = '?typ(Teo2;)^/^^"'"^\ which 
is however different from the average scaling form in Eq. 
( |1 . 1[ ) . The typical autocorrelations now are in a stretched 
exponential form: 



G(r) ^ exp 



1/(^+1) 



6^0, (B7) 



and the relevant scaling variable is given in Eq. (5.10) 
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